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ABSTRACT 

Circular ribbon flares are usually related to spine-fan type magnetic topology 
containing null-points. In this paper, we investigate an X-class circular ribbon 
flare on 2012 October 23, using the multi-wavelength data from the Solar Dy¬ 
namics Observatory, Hinode, and the Ramaty High Energy Solar Spectroscopic 
Imager. In Ca II H emission, the flare showed three ribbons with two highly elon¬ 
gated ones inside and outside a quasi-circular one, respectively. A hot channel 
was displayed in the extreme ultraviolet (EUV) emissions that infers the existence 
of a magnetic flux rope. Two hard X-ray (HXR) sources in the 12-25 keV energy 
band were located at the footpoints of this hot channel. Using a nonlinear force- 
free magnetic field extrapolation, we identify three topological structures: (1) a 
3D null-point, (2) a flux rope below the fan of the null-point, and (3) a large-scale 
quasi-separatrix layers (QSL) induced by the quadrupolar-like magnetic field of 
the active region. We find that the null-point is embedded within the large-scale 
QSL. In our case, all three identified topological structures must be considered 
to explain all the emission features associated with the observed flare. Besides, 
the HXR sources are regarded as the consequence of the reconnection within or 
near the border of the flux rope. 

Subject headings: Sun: flares — Sun: magnetic fields — Sun: UV radiation — 

Sun: X-rays, gamma rays 


1. Introduction 

Flares are believed to be caused by magnetic reconnection that releases part of the 
energy stored in a non-potential field (e.g., see review by Shibata & Magara 2011 and refer¬ 
ences therein). The morphology and evolution of flare ribbons reflect the process of magnetic 
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reconnection. In the classical two-dimensional (2D) flare model (CSHKP; Carmichael 1964; 
Sturrock 1968; Hirayama 1974; Kopp & Pneuman 1976), the flare ribbons appear as a conse¬ 
quence of heating of the lower dense atmosphere caused by the impact of energetic particles 
that are accelerated by magnetic reconnection and propagate downward along the field lines. 
As magnetic reconnection going on, new flare loops are formed with their altitude moving 
up and their footpoints separating from each other perpendicularly to the magnetic polar¬ 
ity inversion line. However, all flares are actually three-dimensional (3D) in nature (e.g., 
Janvier et al. 2013, 2014). For example, magnetic configurations containing 3D null-points 
are observed in many events (Luoni et al. 2007; Masson et al. 2009; Sun et al. 2012a, 2013; 
Mandrini et al. 2006, 2014) in which the flare ribbons, produced by the 3D null-point recon¬ 
nection, are differently shaped from the 2D scenario with two straight ribbons elongating in 
the third dimension. 

Generally speaking, a null-point configuration is embedded in a multipolar magnetic 
field and induces a spine-fan structure around it (Priest & Titov 1996; Moreno-Insertis et al. 
2008; Sun et al. 2012a, 2013; Mandrini et al. 2014). The behavior of the magnetic field close 
to the null-point can be described by the Jacobian matrix of the field. The eigenvalues and 
the eigenvectors of the Jacobian matrix are used to determine the skeleton of the null-point: 
fan plane and spine. For a Jacobian matrix having three eigenvectors with three real and 
distinct eigenvalues (Parnell et al. 1996; Parnell et al. 1997) 3 , the three eigenvalues have 
different signs, one being positive and the other two being negative or on the contrary. This 
is the consequence of the divergence-free condition, which imposes that the sum of these 
three eigenvalues vanishes. The two field lines extending from the null-point in the direction 
of the singly signed eigenvector are defined as the spines, while the field lines extending on 
a plane and corresponding to the other two eigenvectors with the same sign constitute the 
fan separatrix, respectively. The fan surface divides the field into two domains with an inner 
spine and an outer spine. The field-line linkage is discontinuous at the fan separatrix, where 
intense current is prone to be induced (Antiochos et al. 2002; Titov et al. 2003; Galsgaard 
et al. 2003; Pontin & Galsgaard 2007; Pontin & Huang 2012). 

Magnetic reconnection can also occur at other topological structures such as separators 
(Parnell et al. 2010a,b), bald-patches (Titov et al. 1993; Aulanier et al. 1998; Delannee & 
Aulanier 1999; Pariat et al. 2004), or more generally, at quasi-separatrix layers (see Milano 
et al. (1999); Aulanier et al. (2005); Wilmot-Smith et al. (2010); Janvier et al. (2013) for 
MHD simulations, and Lawrence & Gekelman (2009); Gekelman et al. (2012) for laboratory 
experiments). Quasi-separatrix layers (QSLs) have been successfully associated with both 


3 In the other two cases, the eigenvalues are either complex or repeated. Both of them can also determine 
the skeleton structure (Parnell et al. 1996; Parnell et al. 1997). 
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confined and eruptive flares (Demoulin et al. 1997; Savcheva et al. 2012; Guo et al. 2013). 
They are the 3D generalization of separatrices and correspond to volumes of strong gradients 
of the magnetic connectivity (Demoulin et al. 1996a,b; Mandrini et al. 1997; Titov et al. 2002; 
Demoulin 2006). Therefore, the field line connectivities in the two sides of the QSLs are very 
different. Usually, the QSLs are narrow layers. In particular, the thickness of QSLs in 
a flaring configuration is expected to be of the order of 1 km (Demoulin 2006), which is 
beyond the resolution limit in both observations and computations. For viewing purpose, 
one can select a lower threshold of the squashing factor Q, which is a metric of the gradient 
of the magnetic field line mapping (Titov et al. 2002), that yields a finite thickness of QSLs. 
For example, in some cases containing a null-point structure, the spine and fan footpoints 
are shown to be surrounded by QSLs of a finite thickness (Masson et al. 2009, 2012; Sun 
et al. 2013). In observations, flare ribbons always correspond to the intersection of QSLs 
with the upper photosphere (Mandrini et al. 1991; Demoulin et al. 1994). Thus, tracing 
the morphology and evolution of the flare ribbons provides clues to the reconnection process 
(Gorbachev et al. 1988). 

In observations, null-point-related flares would show one or several of the following 
features. First, circular ribbons that are associated with 3D null-point would appear at the 
intersection of the fan with the photoshpere (Liu et al. 2006; Masson et al. 2009; Reid et al. 
2012; Wang & Liu 2012; Sun et al. 2013). Second, if both the inner and the outer spines 
of the null are rooted to the photosphere, a short central ribbon and a remote brightening, 
corresponding to the footprints of the spines, would appear inside and outside the circular 
ribbon, respectively (Liu et al. 2006; Masson et al. 2009; Reid et al. 2012; Wang & Liu 2012; 
Sun et al. 2013). Two other well-known types of eruptive events could also be observed if 
a twisted flux rope is present under the fan of the null. When the flux rope under the fan 
structure loses its equilibrium, its eruption can generate a blowout jet (Moore et al. 2010; 
Sterling et al. 2010; Sun et al. 2012a; Schmieder et al. 2013). This would happen when the 
reconnection in the structure associated with the null-point leads to the disruption of the flux 
rope (Sun et al. 2012a; Schmieder et al. 2013). In such a case, one would observe untwisting 
plasma motions as the flux rope’s disruption results in the launch of non-linear torsional 
Aflven waves propagating along the newly reconnected field lines (Schmieder et al. 2013; 
Pariat et al. 2009, 2010, 2015). On the other hand, if the flux rope maintains its integrity, 
then the null-point reconnection may drive a coronal mass ejection providing that the flux 
rope successfully erupts (Antiochos et al. 1999; Deng et al. 2005; Lynch et al. 2008). 

Recently, the magnetic flux rope has been found likely to appear as a hot channel struc¬ 
ture (Ding et al. 2012; Cheng et al. 2012, 2013, 2014a,b). The hot channel structure and 
the prominence correspond to the main body and the lower part of the flux rope, respec¬ 
tively. Caused by the quick expansion of the flux rope, the upper part of the hot channel 
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seems to separate from the associated prominence. Meanwhile, both of the hot channel and 
prominence experience the same morphology transformation (Cheng et al. 2014a). A 3D 
magnetohydrodynamic simulation of a coronal flux rope by Fan (2012) also supports the 
above observations. In the simulation, the 3D morphology of the hot channel is outlined by 
the isosurface of the temperature at 1.2 MK (Figure 7 in Fan 2012). The hot channel is on 
the top of the central vertical current layer. Thus, the hot channel structure is a signature 
that can be used to trace erupting flux ropes even when there is no filament/prominence, 
which requires further material condensation (Xia et al. 2011, 2012, 2014a,b). 

In this paper, we study the emission features and, in particular, the magnetic topology 
of a circular ribbon flare that occurred on 2012 October 23. We use a nonlinear force-free 
field extrapolation technique to get the 3D magnetic field over the flaring region and further 
deduce the 3D squashing factor Q with the state-of-the-art method proposed in Pariat & 
Demoulin (2012). The magnetic field does contain a null-point and a flux rope below the fan 
structure. We use the magnetic topology to interpret the multi-wavelength emission features 
of the flare, especially its flare ribbons and hard X-ray sources. This paper is organized as 
follows. The observations at different wavelengths are introduced and discussed in Section 2. 
The magnetic field observations and its topology related to the flare are described in Section 
3, followed by a discussion and conclusion in Section 4. 


2. Multi-wavelength Observations of the Flare 
2.1. Call H Observations 

On 2012 October 23, an XI.8 class flare occurred in the active region 11598 (S10, E42) 
that contained two sunspots. In GOES observations, the flare started at 03:13:00 UT, peaked 
at 03:17:00 UT, and ended at 03:21:00 UT. For this event, we study the images of Ca II H 
line that come from the Solar Optical Telescope (SOT; Tsuneta et al. 2008) on board the 
Hinode satellite. The data have a spatial scale of ~ 0.1" per pixel. Since the coordinate 
systems of the SOT and other instruments are not consistent with each other, we align the 
Ca II H image with a shift of (25.9", 20.4") by calculating the cross-correlation between the 
SOT Ca II H image and AIA 1600 A image. 

The central ribbon (Rl) and the western part of the quasi-circular ribbon (R2) bright¬ 
ened first at about 03:14:40 UT (Figure 1). At the initial time, the two ribbons were very 
close to each other so that they were visually seen as one ribbon (Figures la-b). However, 
after carefully checking the two ribbons in the Ca II H movie with high temporal resolution 
(attached to Figure 1), we find that the two ribbons were not connected, but just separated 
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at a place (marked by the circular symbol in Figure la), with the northern part belonging to 
the circular ribbon and the southern part to the central ribbon. The two ribbons, with their 
footpoints close to the polarity inversion line, suggest that the initial magnetic reconnection 
occurred at a relatively low level of the atmosphere. With time going on, the ribbons sepa¬ 
rated from each other with the reconnection site moving up gradually (Figures lb-d and the 
movie of Ca II H), just like what is predicted in the CSHKP model. After the appearance 
of both the central and quasi-circular ribbons, the remote brightening (R3, Figures lc-d) 
started to brighten, with a time delay of less than 20 s that is the time resolution of the SOT 
Ca II H observations. The SOT Ca II H images showed three ribbons around the peak time 
of the flare (Figures lc). The observations displayed some typical features of the circular 
ribbon flare. There existed three brightening parts, a central ribbon (Rl), a quasi-circular 
ribbon (R2) curved around the former, and a remote brightening (R3, c.f., Masson et al. 
2009; Wang & Liu 2012; Deng et al. 2013; Sun et al. 2013; Liu et al. 2013). However, in our 
event, the remote brightening was quite elongated and seemed to be a third ribbon. 

We extract the lightcurves of the Ca II H ribbons by averaging the intensity over the 
ribbons and normalizing them by the pre-flare values, with the saturated values being cor¬ 
rected using a linear interpolation. The selected regions for each ribbon are based on the 
20% contour level of the peak value at 03:22:31 UT. In particular, when the two ribbons 
are close to each other, we use the magnetic polarities to discriminate them and ensure that 
the selected regions do not extend out each ribbon. To estimate the uncertainties of the 
lightcurves , we select different regions ten times and calculate the standard deviations that 
are found to be less than 6% of the mean values. From the lightcurves of the flare ribbons 
(Figure 2), the integrated emissions of the three ribbons display a time sequence in the rise 
phase: the first is the central ribbon, followed by the quasi-circular ribbon, and the last is 
the remote ribbon. Though the emission of the remote ribbon rises last, it gets to peak 
earlier than the other two ribbons. Such a brightening sequence is due to the process of the 
flare reconnection as will be discussed in Section 4. 


2.2. Extreme Ultraviolet Observations 

The 304 A, 171 A and 94 A images from the Atmospheric Imaging Assembly (AIA; 
Lemen et al. 2012) on board the Solar Dynamics Observatory ( SDO ) give a different view 
of this flare from the Ca II H images. The images of the AIA 304 A passband display 
emission features from a low characteristic temperature of 5.0 x 10 4 K, which includes the 
He II line formed in the chromosphere and the transition region. The Fe IX line emission is 
contained in the AIA 171 A passband that originates from the quiet corona and the upper 
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transition region, with a characteristic temperature of 6.3 x 10 5 K. The AIA 94 A passband, 
containing the Fe XVIII line emission, originates from the flaring region corona, with a 
characteristic temperature of 6.3 x 10 6 K. All the observations in these passbands (Figure 
3) with characteristic temperatures ranging over 2 orders of magnitude show the evolution 
of the flare in the higher atmosphere than the Ca II H images. 

One can find that at the initial time, a bright arcade-like structure (ALS) lay low with its 
two footpoints anchored in the solar surface (Figures 3a, 3e, and 3i). From the observations 
of AIA 94 A, three other structures can be found there: highly sheared field lines just above 
the ALS, a dome-like structure over the ALS and the sheared field lines, a loop-like structure 
that connected the dome-like structure and the western part of the solar surface. After a 
short time, at about 03:15:00 UT, this bright ALS rose suddenly and got much brighter 
(Figures 3b, 3f, 3j, and the movie attached to Figure 3). Meanwhile, the sheared field lines 
were squeezed by the rising ALS, and their northern ends also got brighter. Moreover, the 
rising ALS seems to have a little twist, as shown especially in Figure 3f and the movie (the 
AIA 171 A emission) attached to Figure 3, implying that it may be a flux rope. With time 
going on, at about 03:15:25 UT, the rising ALS got higher and it then pushed the dome-like 
structure over it, as revealed in the movie (the AIA 94 A emission) attached to Figure 3. 
Simultaneously with this process, the circular ribbon appeared, which, however, can only 
be seen in the high temperature emissions initially (Figure 3k). Such an evolution process 
suggests that the squeezing of the dome-like structure by the ALS may cause a magnetic 
reconnection between them, and then the accelerated electrons propagate along the field 
lines of the dome-like structure to heat the lower atmosphere that forms the circular ribbon. 
As the ALS rose further up, it became less bright, while the remote and the circular ribbons 
can be seen in the images of all the three wavebands (Figures 3d, 3h, and 31). It is worth 
noting that during the process of the ALS rising, the bright loop-like structure was dragged 
to a higher place (orange dotted line in Figures 3i—1). Then, the remote ribbon appeared 
just at the western end of this loop-like structure at about 03:15:50 UT (Figure 31). This 
implies that energetic electrons are transported through the field lines in this large loop-like 
structure to the lower atmosphere and the remote ribbon is caused by the collision of the 
energetic electrons with the chromospheric plasma. In addition, with the ALS rising, its 
southern footpoint experienced a slight movement but was still anchored inside the circular 
ribbon. There are two possibilities for this motion: one is the lateral expansion of the flux 
rope during its eruption; the other is magnetic reconnection of the flux rope with the ambient 
field that causes an apparent motion of the footpoints. The evolution of the flare in AIA 94 
A is also shown by the movie attached to Figure 3. Comparing the highly sheared loops in 
AIA 94 A at the initial time and the flare loops in Ca II H at a later time, one can find that 
the shear becomes weaker after the flare (Figure 1 and the Ca II H movie attached to it). 
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Such a weakening of the shear angle in a flare process has been studied by many authors 
(Asai et al. 2003; Liu et al. 2010; Inoue et al. 2011; Warren et al. 2011; Savage et al. 2012; 
Aulanier et al. 2012). 

The bright ALS is identified as a hot channel that likely corresponds to a flux rope (Fan 
2012; Cheng et al. 2014a). More evidence for the existence of flux rope can be found in Section 
3. The initial rising of the ALS might be a consequence of tether-cutting type reconnection 
(Moore et al. 2001), breakout reconnection (Antiochos et al. 1999), torus instability (Kliem 
& Torok 2006), kink instability (Torok & Kliem 2005; Guo et al. 2010), flux emergence 
(Chen & Shibata 2000), or magnetic flux cancellation (Wang & Shi 1993; Burtseva & Petrie 
2013). This dynamic process drove the flare to proceed quickly; the 304 A, 171 A, and 94 A 
emissions reached to the peak in about 2 minutes and caused the saturation of the detectors. 
Moreover, as shown in the AIA 94 A movie, with the ALS rising, it became fainter and 
fainter and finally hard to be seen. Then, it is hard to judge what happened to the flux rope 
in the later stage. It might undergo a full eruption or only part of it erupted. Nevertheless, 
no coronal mass ejection was observed during the flare by the Large Angle and Spectrometric 
Coronagraph C2 on the Solar and Heliospheric Observatory. 


2.3. Hard X-ray Emission Features 

Hard X-ray observations of the flare were recorded by the Ramaty High Energy Solar 
Spectroscopic Imager ( RHESSI] Lin et al. 2002) and started at 03:14:36 UT. We reconstruct 
the hard X-ray sources by the Clean algorithm (Hogbom 1974; Hurford et al. 2002) with IF, 
3F, 4F, 5F, 6F, 7F, and 8F detectors for two energy bands of 12-25 keV and 50-100 keV, 
which are superimposed on the Ca II H, 304 A, 171 A, and 94 A images in Figure 1 and 
Figure 3. Two HXR sources can be seen in the 12-25 keV energy band in some of these 
images. The fluxes of the two HXR sources XI and X2 are calculated by summing all the 
pixels in the white boxes (Figure 2a) and then normalized to their peak values (Figure 2b). 

Contrary to the Ca II H emission that was ribbon-like, the HXR sources were localized 
at some particular sites, as usual (Asai et al. 2002; Temmer et al. 2007; Miklenic et al. 
2007; Guo et al. 2012). The source XI was located near the northern end of the parallel 
part of the central and quasi-circular Ca II H ribbons, while source X2 was in the southern 
part of the central ribbon (Figure lc). The HXR source in the 50-100 keV energy band 
overlay the parallel part of the central and quasi-circular ribbons. What is more interesting 
is that this HXR source showed an obvious southward movement along the elongated parallel 
ribbons (Figures lb-d), which is somewhat different from the scenario found by Reid et al. 
(2012), in which the main HXR source, associated with the inner spine and the fan surface, 



keeps stationary for about 60 s, although the secondary source is actually not stationary. 
Meanwhile, the two HXR sources in 12-25 keV also experienced a movement along the 
central ribbon (Figure 2a). The source X2 moved about 10" toward south, while the motion 
of source XI was less than 5" toward south. Note that the cases of both the HXR sources 
moving parallel to the polarity inversion line are not rare, as revealed in a statistic study by 
Bogachev et al. (2005). The authors explained such cases as being caused by the reconnection 
successively occurring along the separator. In those cases, the distance of the two HXR 
sources do not change obviously with time. However, the case here is somewhat different in 
that the distance of the two HXR sources changed obviously with time, with the northern 
one being nearly static but the southern one moving to south. 

In order to make a proper interpretation for the two HXR sources, we further compare 
the HXR images with the EUV emissions at 304 A, 171 A, and 94 A (Figure 3). The result 
shows that two HXR sources correspond to the two footpoints of the rising flux rope, as 
judged from both the hot channel structure (outlined by the dotted line in Figure 3) and 
the magnetic field extrapolation (shown in Section 3.2). The northern footpoint of the flux 
rope moved little and had a good correlation with the northern HXR source. The southern 
HXR source moved toward south and so did the southern footpoint of the flux rope in the 
EUV images, although their locations were not exactly cospatial. Such a spatial relationship 
and similar motion patterns suggest that the HXR sources may be caused by the magnetic 
reconnection within or near the border of the flux rope. Besides, the fluxes of the two HXR 
sources were not peaked simultaneously. The peak of the source X2 was delayed by about 1 
minute with respect to the peak of source XI (Figure 2b). 


3. Magnetic Field Topology of the Flare 
3.1. Magnetic Field Observation and Extrapolation 

The photospheric vector magnetic field data of the active region were observed by the 
Helioseismic and Magnetic Imager (HMI; Schou et al. 2012; Scherrer et al. 2012) on board 
SDO. The data have a pixel spacing of ~ 0.5". As shown in Figure 4, the magnetic flux 
distribution of the active region is, to a first order, equivalent to typical quadrupolar magnetic 
fields, with the four poles defined as PI, P2, Nl, and N2, respectively. The positive pole PI 
was located on the western side, far from the other three poles on the eastern side (P2 lies 
in between Nl and N2). 

We extrapolate the coronal magnetic field by the nonlinear force-free field model with 
the optimization method (Wheatland et al. 2000; Wiegelmann 2004) using the magnetogram 
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observed at 02:59:58 UT before the flare (Figure 4a). First, the 180° ambiguity of the 
transverse components of the magnetic field is removed with the minimum energy method 
(Metcalf 1994; Metcalf et al. 2006; Leka et al. 2009). Since the active region is not located 
close to the disk center, the projection effect should be corrected with the method proposed 
by Gary & Hagyard (1990). As a consequence of removing the projection effect, the geometry 
of the field of view is changed (Figure 4b). We recut the edges to get a rectangle boundary 
for the extrapolation. The extrapolated area is resolved by 200 x 178 grid points with 
8x = 8y ~ 1". Moreover, considering that the force-free and torque-free conditions are 
usually not satisfied for the observed photospheric magnetic field, a preprocessing is applied 
to remove the net force and torque on the boundary (Wiegelmann et al. 2006). The vector 
magnetic field at the bottom boundary after the above preprocessing is shown in Figure 4c. 


Selected field lines from the resulting extrapolation are shown in Figure 5. In order 
to check the force-free state of the extrapolated result, we choose the force-freeness metric 
used by Wheatland et al. (2000), the sine of the angle between the current density and the 
magnetic field averaged by current magnitude: 


(CWsind) = 


Tn Jisinflt 
Si Ji 


where 


sinb, = 


Ji x Eh 

J, B, 


( 1 ) 


( 2 ) 


The metrics of the force-freeness of the extrapolation for the whole 3D computational domain 
and the box containing the flux rope (shown in Figure 5c) are ~ 0.41 and ~ 0.19, corre¬ 
sponding to a current weighted angle of ~ 24° and ~ 11°, respectively. The departure of the 
result from the force-freeness might be caused by the boundary condition, which does not 
satisfy the force-free condition strictly, even after the preprocessing. Adoption of a planar 
geometry of the bottom boundary may also have some influence, since the computational 
domain is quite large and the curvature effect cannot be fully neglected. Nevertheless, the 
current weighted angles in our extrapolation are similar to previous studies (e.g., Schrijver 
et al. 2008; Guo et al. 2010; Sun et al. 2012b). 


3.2. Topological Analysis 

Using the extrapolated 3D coronal magnetic field, we solve the equation Bi(x, y, z) = 0 
(where i = x,y,z) with the modified Powell hybrid method 4 . The magnetic field in the 


4 http: //www.lesia.obspm.fr/fromage/ 



- 10 - 


vicinity of the null is described as the linear term B = M • r, where M is the Jacobian 
matrix for the magnetic field (= dBi/dxj ), and r is the vector of the position. Only 
one null-point is found in the domain, which is marked by a white star in Figure 5. The 
skeleton structure of the null-point, to say, the field lines integrated along the eigenvectors of 
the Jacobian matrix from the null-point, is shown as the red field lines in Figure 5. Next, we 
calculate the squashing factor Q with the method proposed by Pariat & Demoulin (2012). 
This parameter is a geometric measure of the deformation of the elementary flux tube cross 
section and gives the most important information on the magnetic field connectivity gradient. 
The region with high but not infinite Q values is usually a 3D volume with finite thickness. 
By comparison, the separatrices belong to a limit case with the squashing factor Q going 
to infinity and the thickness being infinitesimal. The QSL contains both the regions with 
drastic (2 <C Q ^ oo) and discontinuous (Q = oo) connectivities of magnetic field lines. We 
also calculated the QSLs derived from a potential field extrapolation to emphasize some of 
the identified QSL structures and their interconnections (Figure 5f and Figure 6). 

The topology of the magnetic field reveals three structures as can be seen in Figure 5 
and Figure 6. First, a large-scale QSL, induced by the large quadrupolar-like magnetic fields, 
appears like the cases studied by Demoulin et al. (1996a) and Aulanier et al. (2005), in which 
the footprints of the QSL have two thin strips over two magnetic polarities, respectively. 
Second, we find a single 3D null-point embedded within the large-scale QSL. Third, a flux 
rope is present under the fan structure. The footprints of the large-scale QSL are two C- 
shaped and elongated regions (Cl and C2), corresponding to the footpoints of the yellow 
field lines (Figures 5a-b) and also appearing in the squashing factor map (Figures 5e-f and 
Figure 6). The single null-point is located above a place to the south of the positive pole 
P2, the inner spine is rooted in the south of P2, and the outer spine is rooted in a parasitic 
positive pole to the south of the negative pole Nl. The fan structure divides the volume into 
inner and outer domains, and the intersection of the fan with the lower boundary forms a 
quasi-circular morphology (fan QSL, Figures 5e-f) that is associated with the quasi-circular 
ribbon in AIA 1600 A. For a better comparison, the AIA 1600 A image has been processed 
by removing the projection effect (Figures 5b and 6) with the method proposed by Gary & 
Hagyard (1990). The footprint of the fan structure divides the Cl structure into two parts, 
an inner one and an outer one (Figures 5e-f), which correspond to the inner ribbon and 
the remote ribbon in 1600 A, respectively. The eastern C-shaped region (C2) in the Q map 
(Figures 5e-f and Figure 6) corresponds to the weak ribbon appearing in AIA 1600 A (Figure 
5b and Figure 6) but hardly seen in Ca II H observations (Figure 1). The flux rope lies along 
the western part of the P2-N1 polarity inversion line. Moreover, the two hard X-ray sources 
in 12-25 keV are located beneath the flux rope (Figure 5c). 

To view more clearly the QSL topology, we plot in Figure 7 the distribution of the 
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squashing factor Q in the x — y plane for 2 = 0 and in the x — z plane for different values 
of y. It is seen that viewed in the x — z plane, the squashing factor Q shows a fan structure 
(Figures 7b-c) and also a second dome-like structure above the fan of the null-point (denoted 
as L2 in Figures 7h). Compared with the 3D view of the Q halo (Figure 6 and the movie 
attached to Figure 6), L2 corresponds to the cut in the x — z plane of eastern structure 
(S2) of the QSL induced by the large-scale quadrupolar-like magnetic field. L2 is partially 
coincident with the fan as shown in Figure 7b-d. On the other hand, a vertical structure 
in the x — z plane (denoted as LI in Figure 7h) corresponds to the western structure (SI) 
of the QSL. And the inner domain of the fan consists of two parts, the eastern part (EP) 
and the western part (WP). The flux rope lies in the latter. Such a structure of double 
parts is similar to the ‘fish-bone-like’ structure induced by a magnetic null line (Wang et al. 
2014); however, we do not find a null-line in the magnetic structure of this event. On the 
other hand, we can see that the single null-point is embedded in the western structure of the 
large-scale QSL (SI in Figure 6 and LI in Figure 7i). Thus, the basic magnetic topology of 
this active region is a null-point with its related skeleton embedded in a large-scale QSL and 
a flux rope lying below the null-point-related fan structure. 

Since there is a possibility that the flux rope finally erupts, we check the magnetic 
structure after the flare with NLFFF extrapolation and find that the twisted flux rope seems 
to be replaced by a sheared magnetic field (Figure 8a). The 3D null-point also disappears. 
However, the double C-shaped footprint of the large-scale QSL still remains (Cl and C2), 
as shown in both the Q maps of the NLFFF and potential field (Figures 8b-c). 


4. Discussion and Conclusion 

The flare studied in this paper was associated with a complex magnetic topology com¬ 
prising a flux rope, a 3D null-point, and a large-scale QSL. The flux rope was initially located 
below the fan of the null-point while the null-point was itself embedded in the large-scale 
QSL. The emission features of the flare in Ca II H line, 304 A, 171 A, and 94 A, and HXR are 
also presented in order to study its occurrence and evolution, with regard to the aforemen¬ 
tioned complex magnetic topology. Here, we highlight three points: first, we speculate that 
the HXR sources are caused by the reconnection within or near the border of the flux rope; 
secondly, the evolution of the inner and circular ribbons is associated successively with the 
footprints of the flux rope QSL and the spine-fan QSL; thirdly, a dome-like and a loop-like 
structures in AIA 94 A are observed to correspond to the 3D large-scale QSL halo. 

We identify the ALS (Figure 3) as the flux rope by comparing it with the extrapolated 
magnetic field (Figure 5). Further evidence for the flux rope identification is the hot channel 
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structure in the 304 A, 171 A, and 94 A images (Figure 3). The hot channel structure can 
be a part of the magnetic flux rope (Fan 2012; Ding et al. 2012; Cheng et al. 2012, 2013, 
2014a,b). 

In order to give a proper interpretation for the HXR emission, we make a comparison 
of the HXR sources with the 304 A, 171 A, and 94 A images. The result shows that the 
two HXR sources were located near the footpoints of the flux rope. Further analysis of 
the extrapolated 3D magnetic field (Figures 5c) confirms the spatial relationship between 
the HXR sources and the flux rope footpoints. In particular, with the flux rope rising, 
its footpoints experienced a slight movement (Figure 3) and the HXR sources showed a 
southward movement along the polarity inversion line (Figure 2). The possible origins of 
this slight movement are discussed in Section 2.3. Regardless of these motions, the two HXR 
sources were still located near the footpoints of the flux rope. This is to say, although they 
are not exactly cospatial, their motion patterns are similar. Therefore, the HXR sources are 
considered to be related to the flux rope footpoints. 

A number of previous studies have been done to investigate the magnetic reconnection 
of a 3D null-point, for example, ideal kinematic models (Klapper et al. 1996; Bulanov & 
Sakai 1997; Mellor et al. 2003; Pontin 2011), torsional spine and fan reconnection model 
that occurs due to a rotational disturbance at the fan plane or around the spine (Rickard 
& Titov 1996; Pontin & Galsgaard 2007; Pontin 2011), and the most common spine-fan 
reconnection when a shear disturbance of either the spine or the fan occurs (Pontin et al. 
2007; Pontin 2011). In the above models, the magnetic null point collapses and generates 
electric current locally. The energetic electrons are likely to flow from the null-point along 
the inner/outer spines and the eigenvector with the largest absolute eigenvalue in the fan 
of the magnetic field. However, magnetic reconnection does not solely occur at the null- 
point, but also possibly in extended places with large Q values, like slipping or slip-running 
reconnection (Aulanier et al. 2006; Masson et al. 2009; Reid et al. 2012; Janvier et al. 2013). 
Comparing the skeleton structure of the null-point and the QSL structure with the Ca II H 
emission for the present event, we find that both the footpoints of the inner and the outer 
spines do not correspond to the major bright region; they he even far from, for example, 
the 1600 A emission (Figure 6) and, in particular, the HXR emission (Figure 5c). Figure 6 
shows that the inner spine is rooted in the southern end of the central ribbon. And the outer 
spine is rooted in the eastern end of the extremely elongated remote ribbon. Nevertheless, 
the QSLs correspond to the ribbons well. 

In addition, we find that not only the ribbons in AIA 1600 A correspond to the inter¬ 
sections of the QSLs with the lower atmosphere, but also those in AIA 94 A correspond 
to the QSL halo. By comparing the 3D QSL configuration (Figure 6) with the AIA 94 A 



emission (Figures 3i-l), we find that the field lines embedded in the structure S2 correspond 
to the dome-like struture and those embedded in the structure Si correspond to the loop¬ 
like structure that connects the dome-like structure and the remote ribbon. Physically, the 
QSLs are places where intense currents may appear (Titov et al. 2003; Galsgaard et al. 2003) 
and thus dissipation can take place (Aulanier et al. 2005; Wilmot-Smith et al. 2010; Janvier 
et al. 2013). Therefore, the plasmas in the QSLs could be heated by the dissipation and 
then appear as a bright structure that can delineate the QSL topology. This can explain 
the observational findings in our event. This event is among the very few revealing both a 
dome-like structure and a loop-like structure in the EUV emission, which is consistent with 
a QSL configuration. 

After the above comprehensive analysis, we come up with a plausible interpretation for 
the emission features and time evolution of this flare, which involves three steps. First, a 
flux rope already lies under the fan structure before the onset of the flare; then, magnetic 
reco nn ection occurs inside or near the border of the large flux rope, although we cannot 
determine where it takes place accurately. The rising flux rope, appearing as a large and 
bright ALS (Figure 3), may further incur magnetic reconnection with the surrounding arcade 
field lines. This process feeds new flux to the flux rope and makes it growing gradually. 
Meanwhile, the electrons are accelerated and spiral downward along the field lines to the feet 
of the flux rope, producing two HXR footpoint sources through thick target bremsstrahlung, 
similarly to the event studied by Guo et al. (2012). However, we cannot make it clear what 
causes the time difference of the two HXR sources, which may probably be caused by the 
variation in magnetic and thermodynamic structures along the flux rope. With the flux 
rope expansion and the reconnection going on, the flux rope footpoints likely experience a 
slight movement parallel to the polarity inversion line that may account for the southward 
movement of the HXR source in the 12-25 keV band (Figure 2). Second, the flux rope 
rising may induce the formation of a current sheet just below it where the reconnection 
results in the two flare ribbons in Ca II H, as is predicted in a 3D extension of the CSHKP 
model proposed by Aulanier et al. (2012). At the initial time, the two ribbons are very close 
to each other appearing just like one ribbon, but they can be distinguished in the movie 
attached to Figure 1. Then, they separate from each other and the western one evolves 
to a quasi-circular one. Their locations trace the inner QSL and the western part of the 
fan, respectively. Therefore, in the lightcurves, the central ribbon rises first while the quasi¬ 
circular ribbon starts to rise about half a minute later. Third, the rising flux rope pushes the 
fan structure and the large-scale QSL above and drives the reconnection between the field 
lines embedded in the QSL and fan. This process is supported by the time sequence of AIA 
94 A observations (Figures 3j—1), which shows a dome-like structure being squeezed by the 
rising flux rope. As a consequence of reconnection, electrons are accelerated and transported 



downward along the field lines in the QSLs to heat the chromospheric plasma. Then, the 
Ca II H emission displays a morphology similar to the QSLs at the bottom boundary. It is 
seen that the rising flux rope not only squeezes the structure S2 of the large-scale QSL that 
corresponds to the dome-like structure but also drags the structure SI of the large-scale QSL 
that corresponds to the loop-like structure in AIA 94 A (Figure 31). Such a dragging process 
induces magnetic reconnection in the large-scale QSL and thus the remote ribbon appears 
to be located just at the western footprints of this QSL. The remote ribbon is thought to be 
heated by energetic electrons transported through the large-scale field lines in the large-scale 
QSL; therefore, it appears some time later than the two ribbons mentioned above. After the 
flare, the magnetic field structure contains no null-point, which suggests that the magnetic 
reconnection destroys it during the flare. 

Nevertheless, we are not sure which mechanisms work here to trigger this event. In view 
of the magnetic reconnection process, this event shows some features consistent partly with 
the tether-cutting model (Moore et al. 2001) and partly with the breakout model (Antiochos 
et al. 1999), but neither can fully explain the event. For simplicity, we can define the essence 
of the tether-cutting model as the core field reconnection (the reconnection inside or near 
the border of the flux rope, e.g., the reconnection in a possible hyperbolic flux tube under 
the flux rope) and that of the breakout model as the null-point reconnection (including the 
QSL reconnection) above the flux rope. Usually, both the two models are possible to work 
in one particular event. In our case, the analysis indicates a joint process that involves the 
flux rope rising, core field reconnection (as in the tether-cutting model), and the null-point 
reconnection above the flux rope (as in the breakout model). However, we cannot figure out 
whether the flux rope rising triggers the reconnection or on the contrary. What we propose is 
that the rising of the flux rope and the reconnection process could give a positive feedback to 
each other. While the magnetic reconnection facilitates the rising by increasing the magnetic 
flux of the flux rope (core field reconnection) and decreasing the magnetic tension force above 
the flux rope (null-point and QSL reconnection), the flux rope rising incurs an environment 
to speed up the magnetic reconnection. 

Note that in this event, the extrapolated field lines are not cospatial with the ribbons 
perfectly (Figure 5b). This can be caused by many reasons such as the projection effect, 
systematic errors of the extrapolation, and the non-force-freeness of photospheric field. More¬ 
over, there also exists an uncertainty in the location of the HXR sources owing to the limited 
spatial resolution of the RHESSI reconstructed images. Therefore, at present we cannot give 
more details related to the flare process. Observations with higher spatial and temporal 
resolution would be required in the future in order to make a better understanding of such 
events containing 3D null-point reconnection. 
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Fig. 1.— (a)-(d) SOT/Hinode Ca II H images showing three flare ribbons: a central ribbon 
(Rl), a quasi-circular ribbon (R2), and a remote ribbon (R3). The red circular symbol and 
the arrow in panel (a) indicate the place where the two ribbons separate from each other. 
HXR sources in the 12-25 keV and 50-100 keV energy bands at the closest time to the Ca II 
H images are overplotted with contour levels of 50% and 80% of the peak flux. A Ca II H 
movie is attached to this figure. 
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Fig. 2.— (a) SOT/Hinode Ca II H image at 03:15:57 UT on 23 October 2012 overlaid with 
the HXR sources at different times in the energy band of 12-25 keV. The solid lines indicate 
the contour levels of 50% and 80% of the HXR peak flux. Blue, green, and red colors trace 
the time sequence of 03:15:36 to 03:16:08 UT, 03:16:24 to 03:16:56 UT, and 03:21:48 to 
03:22:20 UT, respectively. The two white boxes delineate the regions over which the HXR 
light curve are extracted, (b) Ca II H lightcurves of the three flare ribbons and the HXR 
fluxes of the two sources in the 12-25 keV energy band. The errorbars stand for the standard 
deviation of the measured fluxes. 
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Fig. 3.— (a)-(d) AIA 304 A images, (e)-(h) AIA 171 A images, and (i)-(l) AIA 94 A images 
showing the evolution of the flare. The blue and yellow dotted curves in (a)-(l) outline the 
rising arcade-like structure. In panel (1), the yellow curve outlines the visible part in 94 A, 
and the blue one is plotted by referring to the structure in 304 A in panel (d). The orange 
dotted lines in panels (i)-(l) outline the loop-like structure (labelled as LS). The black dotted 
line in panel (i) outlines the dome-like structure (labelled as DS). The arcade-like structure 
is labelled as ALS in panel (k). HXR sources in the 12-25 keV and 50-100 keV energy 
bands are overplotted with contour levels of 50% and 90% of the peak flux. A movie of 94 
A emission is attached to this figure. 
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Fig. 4.— (a) Photospheric line-of-sight magnetic field of the active region 11598 at 02:59:58 
UT observed by HMI /SDO before any preprocessing, (b) The vertical magnetic field after 
removal of the 180° ambiguity in the transverse components of the magnetic field and correc¬ 
tion for the projection effect, (c) The vector field at the bottom boundary after preprocessing 
used for the extrapolation. The gray-scale image stands for the vertical component of the 
magnetic field, while the arrows in panel (c) indicate the magnitude and direction of the 
horizontal component. The two positive poles and two negative poles are defined as PI, P2, 
Nl, and N2, respectively. 















27 



Fig. 5.— (a) Magnetic field fines in the computation box, which depicts the null-point 
(red), the QSLs (yellow), and the flux rope (pink-green-blue). The background shows the 
photospheric vertical magnetic field. The white box shows the computational domain of the 
squashing factor Q. (b) The same as panel (a) with the bottom boundary being replaced 
with the AIA 1600 A image at 03:34:40 UT. (c) A side view of panel (a) without the yellow 
magnetic field fines. HXR sources in the 12-25 keV band, with the time interval of 03:16:00- 
03:16:30 UT and contour levels of 50%, 70%, and 90% of the peak flux, are overlaid on the 
photospheric vertical magnetic field. The white box in panel (c) is used to check the force- 
freeness metric of the flux rope, (d) A side view of panel (a) with the bottom boundary 
being replaced with the Q map. (e)-(f) Maps of the squashing factor Q on the bottom 
boundary from NLFFF field and potential field, respectively. The two C-shaped footprints 
of the large-scale QSL are labelled as Cl and C2, respectively. 
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Fig. 6.— (a)-(d) 3D Views of QSL calculated from the potential field. The image on the 
bottom boundary shows the AIA 1600 A emission at 03:34:40 UT. The red lines represent 
the skeleton structure of the null-point. The western and the eastern structures of the large- 
scale QSL that correspond to the loop-like and dome-like structures in AIA 94 A are labeled 
as SI and S2, respectively. A movie showing the 3D QSL is attached to this figure. 
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Fig. 7.— (a) Map of the squashing factor Q in the x-y plane at z = 0 calculated from the 
NLFFF extrapolation, (b)-(i) Maps of Q slices in the x-z planes with different y values. The 
red star in panel (i) marks the position of the null-point. LI and L2 correspond to the cuts 
in the x — z plane of the western and the eastern structures (SI and S2) of the large-scale 
QSL induced by the quadrupolar magnetic field. 
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Fig. 8.— (a) Magnetic field after the flare calculated from the NLFFF extrapolation, in 
which the red field lines are highly sheared and the yellow field lines belong to the QSL. 
(b)-(c) Maps of squashing factor Q on the bottom boundary calculated from the NLFFF 
extrapolation and the potential field, respectively. The double C-shaped footprints of the 
large-scale QSLs after the flare are labelled as Cl and C2, respectively. 











